Algebraic interface-based coarsening AMG preconditioner for multi-scale sparse matrices with applications to radiation hydrodynamics computation
نویسندگان
چکیده
Coarsening is one of the crucial components of an algebraic multigrid (AMG) method for iteratively solving sparse linear systems arising from scientific and engineering applications. It largely determines the complexity of the AMG iteration operator. Usually, high operator complexities lead to fast converge of the AMG method, however, they require additional memory and as such do not scale as well in parallel computation. On the other hand, though low operator complexities improve parallel scalability, they often lead to deterioration in convergence. This paper introduces a new type of coarsening strategy, called algebraic interface based coarsening that yields better balance between convergence and complexity for a class of multi-scale sparse matrices. Numerical results for various model-type problems and a radiation hydrodynamics practical application are provided to show the effectiveness of the proposed AMG solver.
منابع مشابه
Reducing Complexity in Parallel Algebraic Multigrid Preconditioners
Algebraic multigrid (AMG) is a very efficient iterative solver and preconditioner for large unstructured linear systems. Traditional coarsening schemes for AMG can, however, lead to computational complexity growth as problem size increases, resulting in increased memory use and execution time, and diminished scalability. Two new parallel AMG coarsening schemes are proposed, that are based on so...
متن کاملA fast algebraic multigrid preconditioned conjugate gradient solver
This work presents a new approach for selecting the coarse grids allowing a faster algebraic multigrid (AMG) preconditioned conjugate gradient solver. This approach is based on an appropriate choice of the parameter a considering the matrix density during the coarsening process which implies in a significant reduction in the matrix dimension at all AMG levels. 2005 Elsevier Inc. All rights rese...
متن کاملOn the Scalability of the Albany/FELIX first-order Stokes Approximation ice Sheet Solver for Large-Scale Simulations of the Greenland and Antarctic ice Sheets
We examine the scalability of the recently developed Albany/FELIX finite-element based code for the first-order Stokes momentum balance equations for ice flow. We focus our analysis on the performance of two possible preconditioners for the iterative solution of the sparse linear systems that arise from the discretization of the governing equations: (1) a preconditioner based on the incomplete ...
متن کاملFast hierarchical solvers for sparse matrices
Sparse linear systems are ubiquitous in various scientific computing applications. Inversion of sparse matrices with standard direct solve schemes are prohibitive for large systems due to their quadratic/cubic complexity. Iterative solvers, on the other hand, demonstrate better scalability. However, they suffer from poor convergence rates when used without a preconditioner. There are many preco...
متن کاملDELFT UNIVERSITY OF TECHNOLOGY REPORT 13-11 Evaluation of Multilevel Sequentially Semiseparable Preconditioners on CFD Benchmark Problems Using IFISS
This paper studies a new preconditioning technique for sparse systems arising from discretized partial differential equations (PDEs) in computational fluid dynamics (CFD), which exploit the multilevel sequentially semiseparable (MSSS) structure of the system matrix. MSSS matrix computations give a data-sparse way to approximate the LU factorization of a sparse matrix from discretized PDEs in li...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Numerical Lin. Alg. with Applic.
دوره 24 شماره
صفحات -
تاریخ انتشار 2017